Fermi-liquid based theory for the in-plane magnetic anisotropy 
in untwinned high-Tc superconductors 
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Using a generalized RPA-type theory we calculate the in-plane anisotropy of the magnetic excita- 
tions in hole-doped high-Tc superconductors. Extending our earlier Fermi-liquid based studies on 
the resonance peak by inclusion of orthorhombicity we still find two-dimensional spin excitations, 
however, being strongly anisotropic. This reflects the underlying anisotropy of the hopping matrix 
elements and of the resultant superconducting gap function. We compare our calculations with new 
experimental data on fully untwinned YBa2Cu3 06.85 and find good agreement. Our results are in 
contrast to earlier interpretations on the in-plane anisotropy in terms of stripes (H. Mook et at., 
Nature 404, 729 (2000)), but reveal a conventional solution to this important problem. 



Since the discovery of high-Tc superconductors, its 
mechanism is still under debate. Perhaps one of the most 
important question concerns the role played by spin ex- 
citations in these materials. For example, one scenario 
of superconductivity in layered cuprates suggests that 
Cooper-pairing is due to an exchange of antiferromag- 
netic spin fluctuations [1]. In this respect, an under- 
standing of the so-called resonance peak observed by in- 
elastic neutron scattering (INS) experiments [2,3] at the 
antiferromagnetic wave vector Qaf and energy ui « Ures 
plays an important role in the phenomenology of high-Tc 
superconductors. Among various explanations over last 
years there are two most probable scenarios for the for- 
mation of the resonance peak. The first one suggests that 
the two-dimensional Cu02 layers are intrinsically unsta- 
ble towards a stripe formation with one-dimensional spin 
and charge order [4]. In this picture, the resonance exci- 
tations can be interpreted in terms of excitation spectra 
in a bond-centered stripe state with long-range magnetic 
order [5,6]. More technically, the resonance peak in the 
stripe-ordered phase corresponds to a saddle point in the 
dispersion of the magnetic excitations. In the other ap- 
proach that is a conventional Fermi-liquid one, the res- 
onance peak arises as a particle-hole excitation (or spin 
density wave collective mode) in a c?3;2_j,2-wave super- 
conductor and is a result of the strong feedback of the 
superconductivity on the dynamical spin susceptibility 
below Tc [7-12]. 

In order to distinguish between both pictures, a 
detailed analysis of untwinned cuprates is necessary. 
Recently, INS study in the fully untwinned high- 
temperature superconductor YBa2Cu3 6.85 reveals two- 
dimensional character of the magnetic fluctuations [13] in 
contrast to the previous conclusions from measurements 
in the partially untwinned samples [14]. Here, motivated 
by recent experiments we analyze the in-plane anisotropy 
of the magnetic excitations in hole-doped high-Tc super- 
conductors within a conventional Fcrmi-liquid and gen- 



eralized RPA-like approach. Extending our earlier stud- 
ies on the resonance peak by the inclusion of a small 
orthorhombicity we still find two-dimensional spin ex- 
citations, however being strongly anisotropic, reflecting 
the underlying anisotropy of the hopping matrix elements 
and of the resultant superconducting gap function. 

In order to describe the phenomenology of the su- 
perconducting cuprates we employ an effective one-band 
Hubbard Hamiltonian for the Cu02-plane 



(1) 



where is a creation operator of an electron with spin 
a on site i, U denotes the on-site Coulomb repulsion, 
and tij is a hopping matrix element in the Cu02-plane. 
Here, we use the six parameter fit of the energy dis- 
persion suggested in Ref. [15] with the following chem- 
ical potential and hopping amplitudes (^, ti, ^5) (the 
units are in eV): (0.1197, -0.5881, 0.1461, 0.0095, -0.1298, 
0.0069). The lattice constants are set to unity. To de- 
scribe the orthorhombic distortions we introduce a pa- 
rameter Sq which leads to an anisotropy in the hop- 
ping integrals along and perpendicular to the chains in 
YBa2Cu3 06.85(YBCO). This one-band approach seems 
to be justified at least in the optimally-doped cuprates 
because upon hole doping into the Cu02-planc antifer- 
romagnetism disappears due to Zhang-Rice singlet for- 
mation and quenching of Cu spins. Further doping then 
increases the carrier mobility and a system of strongly 
correlated quasiparticles occurs [16]. 

Normal state: Before analyzing the superconducting 
state it is instructive to understand how the normal state 
properties and the electronic structure of a Cu02 plane 
are affected by the presence of the orthorhombic distor- 
tions. In Fig.l we show the calculated density of states 
(DoS) and Fermi surface topology as a function of the 
orthorhombic distortions. 
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FIG. 1. (color online) Calculated density of states with and 
without orthorhombicity. For a comparison the corresponding 
changes of the Fermi surface topology are shown in the inset. 
The arrows refer to the two 2kF instabilities as described in 
the text. 

Without orthorhombicity the DoS reveals a pro- 
nounced van-Hove singularity (VHS) being approxi- 
mately 19meV below the Fermi level in good agreement 
with early ARPES experiments on YBCO [17]. Due to 
the vicinity of the VHS to the Fermi level, the effect 
of the orthorhombic distortions is quite strong. First, 
at So ~ 0.03 the singularity is suppressed and shifted 
slightly above the Fermi level. With further increase of 
Sq it splits into two peaks. Similar changes occur for the 
Fermi surface. As a function of the orthorhombicity its 
topology changes and the Fermi surface closes around the 
(— TT, 0) and (tt, 0) points which also leaves an impression 
that the system turns towards a quasi-one-dimensional 
ones. One of the immediate consequence of these changes 
is that the VHS will be present below the Fermi level only 
around (0, ±7r) points. This is also consistent with more 
recent ARPES data on untwinned YBCO [18] where a 
suppression of the ARPES intensity was observed around 
(±7r, 0) due to absence of the VHS. However, in Ref. [18] 
the change in the Fermi surface topology around (±7r, 0) 
was not confirmed. It is important to note that this 
Fermi surface deformation breaks the point-group sym- 
metry and looks similar to what is expected for the case 
of a (ia;2_j,2-wave Pomeranehuk instability due to strong 
electron-electron interactions [19]. 

What happens to the spin response if the electronic 
properties are changed due to orthorhombicity? We cal- 
culate the real part of the bare spin susceptibility in the 
normal state, 



. s ^ /(£k) - /(£k-Hq) 

^ £k+q - Sk + LU + iO+ 
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where /(sk) is the Fermi function. In Fig. 2 we show the 
calculated Rexo(q, = 0) as a function of the trans- 
ferred momentum q. Without orthorhombicity its peak 
structure reflects two 2\ip instabilities of the Fermi sur- 
face (see inset of Fig.l). The first peak one corresponds 
to the quasi-one-dimensional wave vector connecting the 



Fermi surface (FS) around (0, ±7r) points and the other 
one refers to the wave vector connecting the FS along the 
diagonal of the first Brillouin Zone. Note, due to strong 
nesting of the Fermi surface the second structure has the 
form of a plateau. If the orthorhombicity is present the 
first peak is suppressed and shifted towards higher q val- 
ues. This is due to the fact that the Fermi surface closes 
around (±7r, 0) point removing this instability there and 
moving them to higher q values around (0, ±7r) points. 




FIG. 2. (color online) Calculated Re Xo(ci, 0) along the 
path (0, 0) (tt, 0) (tt, tt) — > (0, 0) for various values of 
the orthorhombic distortions. The arrows indicate the peaks 
in Rexo(q, 0) arising from 2kF instabilities as shown in Fig.l. 

On the other hand the diagonal second 2\s.f instability 
remains mainly unchanged and even became more pro- 
nounced, since the plateau around Qa_f = (tt, tt) is sup- 
pressed due to the changes of the Fermi surface topology 
at the parts connected by this wave vector. 

We safely conclude that an orthorhombic distortion 
change strongly the electronic structure of the Cu02- 
planc and yield characteristic changes of the Fermi sur- 
face topology, quasiparticle density of states, and the spin 
susceptibility. We would like to stress that despite the 
changes of the Fermi surface topology indicating the ten- 
dency towards quasi-one-dimensionality, the static spin 
susceptibility remains mainly two-dimensional. 

Superconducting state: In our one-band model, we fur- 
ther assume that the same quasiparticles are participat- 
ing in the formation of antiferromagnetic fluctuations and 
in Cooper-pairing due to these fluctuations. This leads 
to the generalized Eliashberg equations which have been 
derived and discussed in Refs. [20,21]. Although these 
equations cannot describe the metal-insulator transition 
properly, we would like to stress that they allow us to cal- 
culate all properties of the system self- consistently such 
as the elementary excitations, the superconducting or- 
der parameter, and the dynamical spin susceptibility, for 
example. 

In the pure tetragonal case the resulting superconduct- 
ing order parameter has dr^2_y2-wave symmetry. How- 
ever, in presence of orthorhombicity the superconduct- 
ing order parameter changes, since s-wave and d-wave 
symmetries belong now to the same irreducible represen- 
tation of the point group symmetry. The total supercon- 
ducting gap has the form (weak-coupling limit, Z — 1) 
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A(k) = g{5o)As + /((5o)Arf(k) 
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FIG. 3. (color online) Calculated influence of the or- 
thorhombic distortions (s-wave component of the supercon- 
ducting gap and changes in the electronic structure) on the 
resonance peak. Inset: Influence of the admixture of s-wave 
superconducting gap only on the resonance peak formation in 
the tetragonal system. Here, in order to fit the position of the 
resonance at Qaf around 41meV we use C/=0.155eV. Note, 
we further employ the damping F = 2.4meV. 

where, for simplicity, we employ g{So) ~ So, fiSo) = 
1 — ^0 and Ad = Ao(cosfca: — cosfcy)/2, = Aq. For the 
set of parameters described above, we use Aq = 26meV. 
Note that the additional s-wave component leads to an 
in-plane anisotropy of the gap function and thus to dif- 
ferent maximum gap values between (±7r, 0) and (0, ±7r) 
as observed in Ref. [22]. 

The most interesting question is: what happens to 
the dynamical spin susceptibility in the superconduct- 
ing state if an orthorhombic distortion is present. Below 
Tc, within generalized RPA, the imaginary part of the 
dynamical spin susceptibility is given by 

imxo(q,^^) 



Imx(q,t^) = 



1 - [/Rexo(q, + Unmxliq, u;) ' 



(4) 



where xo is the BCS Lindhard response function [15]. 
Without orthorhombicity, the dx2-y2 superconducting 
gap opens rapidly due to a feedback effect on the ele- 
mentary excitations [23] yielding a jump at 2Ao in Im^o 
and the resonance condition [7-9] 



1 - URe xo(q = Qaf,uj 



= 
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is fulfilled. Since Imxo is zero below 2Ao, the reso- 
nance condition (5) reveals a strong delta-like peak in 
Imx which occurs only below Tc. Note, its position is 
mainly determined by the maximum of the o?-wave su- 
perconducting gap Aq and also by the proximity to an 
antiferromagnetic instability described by the character- 
istic energy scale LUsf (roughly the peak in Im x(Qa_f, i-^) 
in the normal state). Then, the resonance peak scales 
with the maximum of the d-wave superconducting gap 
in optimally doped and overdoped compounds. On the 



other hand, in the underdoped cuprates it rather scales 
with LJsf [7-9] due to stronger antiferromagnetic fluc- 
tuations. Thus, one finds for the whole doping range 
t^res/kBTc ~ const [7-9] in good agreement with experi- 
ments [24]. 

In Fig. 3 we analyze the influence of the orthorhombic 
distortions on the resonance peak. One clearly sees that 
the orthorhombicity slightly shifts the resonance peak to- 
wards higher energies and reduces its intensity for in- 
creasing Sq. As already mentioned, the position of the 
resonance peak is determined by the strength of the an- 
tiferromagnetic fluctuations present in the normal state 
and by the d-wave superconducting gap. Both are de- 
creasing due to orthorhombicity. Namely, as one sees 
from Fig. 2, the susceptibility is decreasing around Qaf 
due to a change of the Fermi surface topology. Thus, 
in the superconducting state the resonance is shifted to- 
wards higher energies. On the other hand, the maximum 
of the d— wave gap is decreasing (as we see from Eq.(3)) 
which would shift the resonance energy towards lower 
values. Most importantly, we see that the deformation of 
the electronic structure and Fermi surface topology due 
to orthorhombicity dominate the effect from considering 
solely than the increase of the s-wave component of the 
superconducting gap (see inset). Noticeable changes oc- 
cur only for very large values of Sq when the electronic 
structure is already strongly distorted. Thus we con- 
clude, the observed slight shift of the resonance peak po- 
sition in untwinned [13] occurs due to the strong changes 
in the electronic structure induced by the orthorhombic 
distortion rather than due to additional s-wave compo- 
nent of the superconducting gap. 
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FIG. 4. Calculated dispersion of the resonance peak to- 
wards lower energies for fixed Qy — n without including the 
effect of the orthorhombicity. 

Comparison with experiment: Let us now turn to the 
broader analysis of the dispersion of the resonance ex- 
citations below ujres- In the tetragonal case, going away 
from the points of the Fermi surface connected by the an- 
tiferromagnetic wave vector C^af , we are moving towards 
the diagonal of the BZ where the d^2_y2-wave supercon- 
ducting gap is zero. As a result the resonance energy, 
LUres shifts towards smaller values. This results in the 
parabolic shape of resonance energy dispersion as shown 
in Fig. 4 calculated for fixed qy ~ tt. This parabolic be- 
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havior obtained in our calculations agrees well with the 
experimental findings of Bourges et al. [24] . 

What is happening below the resonance threshold {ui < 
LUres) in fuUy untwinned YBCO for constant energy scans 
as a function of the momenta qx and qyl In Fig. 5 we 
show the calculated projected momentum dependence 
of Imx(q, = 35meV) without (a) and with (b) or- 
thorhombicity. In accordance with ab-initio calculations, 
we have chosen 60 ~ 0.03 [25]. 
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FIG. 5. (color) Calculated normalized two-dimensional in- 
tensity plot for a constant energy of hu = 35 meV, (a) with- 
out, (b) with inclusion of orthorhombicity {So = 0.03). Note, 
for convenience we have chosen the axes as they are in the 
experimental work of Ref. [13]. 

In the tetragonal case one sees that the spin excita- 
tions form a circle around (tt, tt) with four pronounced 
peaks at (tt ± qo,T^) and {tt,tt ± go)- The origin of the 
peaks is clear: away from Qaf we are connecting points 
at the Fermi surface which lie closer to the diagonal of 
the BZ. The superconducting gap tends to zero there 
and thus the position and the intensity of the resonance 
peak are decreasing. However, for the diagonal wave 
vectors (tt ± qo , tt ± qo) it happens faster than for the 
vector (tt ± qo,Tr) or (tt, tt ± qo). Therefore, effectively 
the latter peaks are 'closer' to the resonance condition 
at CIaf = (""j"'); their intensities are higher than those 
for the other wave vectors. This explains the observed 
symmetry of the dominant spin excitations for uj < iOres 
shown in Fig. 5(a). For the orthorhombic case the situ- 
ation is changing. The ring of the excitations becomes 
distorted and, most importantly, there are only two well 
pronounced peaks. The latter is a result of strongly 
changed electronic properties, in particular, the topology 



of the Fermi surface (see inset of Fig.l for comparison). 
Our main result is that, despite there are only two pro- 
nounced peaks, the resonant spin excitations remains ba- 
sically two-dimensional. This is in good agreement with 
recent experiments [13]. Furthermore, this result based 
on a standard Fermi-liquid approach is in contrast to the 
stripe scenario of the resonance peak [5]. 

In summary, we have analyzed the in-plane magnetic 
anisotropy in high-Tc superconductors with orthorhom- 
bic distortions employing a generalized RPA-type theory 
and compared our results with INS data on fully un- 
twinned YBCO. We find that due to changes in the elec- 
tronic structure and the FS topology the resonance peak 
is slightly shifted towards higher energies and that for 
UJ < LOres the dominant spin excitations form a 2D ring- 
like structure around Qaf with four pronounced peaks 
(tetragonal case) . The orthorhombic distortions suppress 
two of these peaks, however, the overall structure of the 
excitations in Im^ remains two-dimensional which agrees 
well with recent experimental data by V. Hinkov et al. 
[13]. Our results provides an alternative picture based 
on a conventional Fermi-liquid theory in contrast to the 
stripe scenario. 
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